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Abstract. We consider the Fuchsian linear differential equation obtained 
(modulo a prime) for x'^^ i the five-particle contribution to the susceptibility of the 
square lattice Ising model. We show that one can understand the factorization of 
the corresponding linear differential operator from calculations using just a single 
prime. A particular linear combination of x^^^ and x'^' can be removed from x'^' 
and the resulting series is annihilated by a high order globally nilpotent linear 
ODE. The corresponding (minimal order) linear differential operator, of order 29, 
splits into factors of small orders. A fifth order linear differential operator occurs 
as the left-most factor of the "depleted" differential operator and it is shown 
to be equivalent to the symmetric fourth power of Le, the linear differential 
operator corresponding to the elliptic integral E. This result generalizes what we 
have found for the lower order terms x^^^ a-nd x^'^h We conjecture that a linear 
differential operator equivalent to a symmetric (n — l)-th power of Le occurs as 
a left-most factor in the minimal order linear differential operators for all x'"' 's. 
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1. Introduction 

Wu, McCoy, Tracy and Barouch [1] have shown that the magnetic susceptibihty of the 
square lattice Ising model can be expressed as an infinite sum of contributions, known 
as n-particle contributions, so that the high-temperature susceptibility is given by 

kT-XH{w) = ^x(^"+i)H = i-(l-s4)^-5]x^2„+i)(^) 



and the low-temperature susceptibility is given by 

kT-XL{w) = ^^('"^H = (1 - 1/^4)4 •^x^2„)(^) (2) 

in terms of the self-dual temperature variable w — ^s/{l -f s^), with s = smh{2J/kT). 

As is now well known [T], the n-particle contributions have an integral 
representation and are given by the {n — l)-dimensional integrals [U [31 HI [5] 
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valid for small w and, elsewhere, by analytic continuation. The variable w corresponds 
to small values of s as well as large values of s. It is worth noting that the series 
expansions for x'"^ in the variable w have integer coefficients. From the first x^"', 
the coefficients for generic n can be inferred [6] 

x'"' - 2" • u;"' ■ (l + An^ -w^ + 2- {An^ + ISn^ + 1) • 

8 
3 

where the coefficients are valid up to nP' for n > 3, w'* for ?i > 5 and for rt > 7 
(in particular it should be noted that x^"^ is an even function only for n even). 

In previous work ^7] we performed massive computer calculations to obtain the 
susceptibility of the square lattice Ising model and the n-particle contributions x^"' • 
These calculations confirmed previously [HI [5] conjectured singularities for the linear 
ODEs of X*'"'' (for n = 5 and 6) and yielded the values of the associated local 
exponents. In addition some light was shed [7] on important physical problems such 

j The Fermionic term G'"' has several representations [5]. 
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as the existence of a natural boundary for the susceptibihty of the square lattice Ising 
model and the subtle resummation of logarithmic contributions from individual x^'^^^s 
resulting in the power-law behaviour of the full susceptibility x- 

As far as the five-particle contribution to the susceptibility is concerned, a long 
series S{w) for x^^^ was generated modulo the prime p,. = 2^^ — 19 from which we 
obtained [7] the corresponding Fuchsian differential equation. This Fuchsian linear 
ODE is of order 33 and we denote by L33 its linear differential operator. The 
calculation of the series is very time consuming and one cannot calculate (given 
presently available computational resources) the many series modulo various primes 
required to reconstruct, through the Chinese remainder procedure, the exact series 
for x''^\ and, from this, deduce the corresponding exact Fuchsian linear ODE. Our 
purpose here is, using the series and the linear ODE obtained modulo a single prime, 
to perform, as far as possible, the factorization of the linear differential operator L33 
and gain as deep an understanding as possible of the various factors occurring in its 
exact factorization (over the rationals). 

In particular, we find that a certain linear combination of x^^^ and x^'^^ can be 
removed from x^^^ and the resulting series is a solution of an order 29 linear ODE. We 
develop methods which enable us to show that the corresponding linear differential 
operator L29 splits into several factors and we present arguments that the order of 
any individual factor does not exceed five. The factor L5 of maximum order occurs as 
the left-most factor of L29- We show that L5 is eauivaleni fffl to the symmetric fourth 
power oi Le, the linear differential operator corresponding to the complete elliptic 
integral E, see ([SS)) . This result generalizes what we have found in [lOl [HI [12] for the 
lower terms x^'^^ and x^''''- We therefore conjecture that a linear differential operator 
Ln, equivalent to the symmetric (n— l)-th power of Le, occurs as the left-most factor 
in the (minimal order) linear differential operators for all the x''"''s. 

2. Deciphering the structure of x''"''^ direct sums, symmetric powers and 
modular forms 

A linear differential operator L can be viewed formally as a non-commutative 
polynomial in w and D^,, where D^, = d/dw is the derivation (or derivative) with 
respect to w. In previous works [Tni[IIl[Il] we have shown that the (minimal order) 
linear differential operators for x'-'^'' and x'''*'' (called respectively L7 and Lio) have a 
"Russian-doU" structure involving the differential operators Li and Nq for x'-^'' and 
X^^\ respectively. More precisely, x^^'' and x''^'' are solutions of the linear ODEs 
corresponding to L7 and Lio, respectively. In terms of linear differential operators 
this means that Li (resp. A'o) right-divides L7 (resp. Liq). Note that throughout 
this paper when we talk about a homogeneous linear differential equation and its 
associated differential operator we will use the terms ODE and differential operator, 
interchangeably. 

One might then conjecture that this structure extends to the linear differential 
operator L33 (for x^^^) and the finear differential operator L7 (for x^'^^)- We note that 
the singularities for the ODEs corresponding to x*-^-* and x''^'' are consistent with this 
assumption, that is all the singularities of L7 also occur in L33. The check of the right 
division between operators can be done simply by generating the series L7 (S{w)) and 
obtaining the corresponding linear ODE. If the order of this latter linear ODE is less 

ft For the notion of equivalence of linear differential operators see |18II19I . 
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than 33, the assumption is verified, i.e. L7 right-divides L33. As reported in [7] this 
procedure leads to the factorization 

L33 = N26 ■ Lt, (9) 

where N26 is a hnear differential operator of order 26. 

A stronger conjecture amounts to saying that the linear differential operator for 
occurs as part of a direct-sum decomposition of the linear differential operator 
for x^"^^^- This conjecture is based on our findings [TUl El [H] that the linear 
combinations 6x'^"+^-' — nx^"-* with n= 1, 2, happen to verify a linear ODE of lower 
order than the one for x*-"^^-*- The conjecture was verified for n = 3 by obtaining [7] 
the minimal order Fuchsian linear ODE for 6 x^^^ — 3x'-^-' , which happens to be of 
order 30, so that 

^33 = Lt ® L30. (10) 

This direct-sum structure, where a seventh order linear differential operator 
combined with an order 30 linear differential operator gives rise to an order 33 linear 
differential operator, leads to the conclusion that there must be a fourth order linear 
differential operator that right divides both L30 and L7. This kind of direct-sum 
structure (jlOp is not seen in L7 or Lio, the linear differential operators for x^^-* and 

Recalling 13J the factorization of L7, the fourth order linear differential operator 
reads (following the notation of Eq. (7) in [T3]) 

B2 • Ti • Li = B2-0i-Ni = Xi-Z2-Ni ^ Li ® (Z2 • TVi) (11) 

and the factorization of L30 in pH)) becomes 

^30 = ^26 -(il ® iZ2-m)). (12) 

In a further step, NickelH has succeeded in showing that the differential operator Li 
(corresponding to x^^'') occurs via a direct sum in L30. This was done by considering 
the series for the combination 

X^'^ ~lx^'^ -a-x^'l (13) 

If a rational value of a can be found such that the resulting linear ODE has an order 
less than 30 (in fact equal to 29), then the direct sum assumption has been validated. 
This happens with a = —1/120 and the final result is that the combination 

is annihilated by an order 29 linear differential operator that we denote L29, leading 
to conclude that 

^30 = il 2.29. (15) 

At this point, guided by our results for the three terms x'""*; = 3, 4, 5, one may 
wonder whether there is a common structure to the corresponding linear differential 
operators. 

Recall that the ODE for x*^"^"* is of order seven and that x^^'' can be written as 
6 



t We are grateful to B. Nickel for this result. 
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where is a solution of a sixth order linear ODE. We thus have the direct sum 
decomposition 

Lr = Li®Le. (17) 

The sixth order operator Lq has a third order linear differential operator Y3 as a 
left-most factor 

Le ^ Y3-Z2-Ni, (18) 

and we have given the solutions of the linear ODE corresponding to Y3 in [T3] . These 
solutions can be written [13j as a homogeneous polynomial of the complete elliptic 
integrals K and E with homogeneous degree two (the order of ¥3 minus one) . 

Next we consider the tenth order linear ODE for x*-'*-'. Recall that x'^"* can be 
written as 

X^'^ = Ix^'^ + (19) 

where ^'-^^ is a solution of an eighth order linear ODE. We thus have the direct sum 
decomposition 

Lio = No(SLs. (20) 

The eighth order operator Lg, has a fourth order linear differential operator M2 as its 
left-most factor: 

Ls - M2 • i4. (21) 

The fourth order linear differential operator L4 factorizes into four first order linear 
differential operators as shown in Eq. (42) of [11]. As mentioned in 13J (and now 
given explicitly in [Appendix A] ) , the linear ODE corresponding to M2 annihilates a 
homogeneous polynomial of K and E of (homogeneous) degree three, i.e. the order of 
M2 minus one. 

Similarly, we have shown that for x*-^' 

^(5) ^ 1^(3) -^x(^) -t-cl>(^), (22) 

where is annihilated by an order 29 linear ODE whose corresponding differential 
operator we denote as L29. 

We conjecture that once the x*^""* are depleted of the contributions x'-""^'"-' of 
lower index (with coefficients a„^2fc, where a„_2 — {n — 2)/6, and the remaining 
coefficients are to be determined numerically) 

X'"^ = «n-2 • X'"~'^ + «n-4 ' X^"-'^ + • • • + (23) 

the depleted series will be solutions of linear ODEs of minimal order q, whose 
corresponding (minimal order) linear differential operators factorize as 

Lq = Ln ■ ig-nj (24) 

and where the linear ODE corresponding to the left-most factor i„ has as a solution 
a homogeneous polynomial of complete elliptic integrals E and K of degree n — 1 (in 
other words L„ is equivalent to the (n — l)-th symmetric power of L^;, annihilating 
E, see below). 

This is what happens for the terms x'^^'' and x^^"*- O^ic of the purposes of this 
paper is to show that this structure also holds for x*-^-* • This amounts to demonstrating 
the occurrence of a fifth order linear differential operator L5 at the left of L29, with 
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L5 being equivalent to the symmetric fourth power of the linear differential operator 
L E corresponding to the complete elliptic integrals E (or K) . 

Before proceeding to show how we achieved this goal, some general remarks are 
in order. For an integral representation of a function its series expansion S{x) around 
the origin (x — 0) is unique. This series can be annihilated by (is a solution of) many 
linear ODEs of order Q and degree D (see Appendix B): 

Among all these linear ODEs there is one of minimal order q and it is unique (its 
corresponding degree will be denoted by Dq). In terms of linear differential operators, 
the minimal order differential operator appears as a right-factor in the non-minimal 
order linear differential operators. The minimal order linear ODE may contain a 
very large number of apparent singularities and can thus only be determined from a 
very large number of series coefficients (generally speaking [q + l)(I?o + 1) terms are 
needed). Other (non-minimal order) linear ODEs, because they carry polynomials 
of smaller degrees, may require fewer series coefficients in order to be obtained. For 
any Q > q, a, linear ODE annihilating S{x) (i.e. LqdISIx)) = 0), can be founcQ 
for D sufficiently large, and if Q is small enough we can choose Q and D such that 
{Q -|- 1)(-D -|- 1) < {q + l)(-Do + !)■ Among the non- minimal linear ODEs there will 
generally be one requiring the minimal number of terms; in a computational sense one 
may view this as the "optimal linear" In the case of x*-^-* this optimal linear 

ODE can be discovered from some 7400 terms while the minimal order linear ODE 
requires almost 49100 terms. So when we consider, for instance, a linear differential 
operator such as L29 (of minimal order 29), we are dealing in fact (as far as the 
computations are concerned) with a much higher order linear differential operator. 

The knowledge about the minimal order is "inferred" from many non-minimal 
order ODEs by using the remarkable formula (|26p below, which we reported in [7]. 
Seeking a Fuchsian linear ODE of order Q and degree D which annihilates a given 
series requires a certain number of coefficients N. Formula (|26p relates this number 
of required coefficients N to the order Q and degree D of the Fuchsian linear ODE. 
Remarkably, this "ODE formula" gives the number of required coefficients as a 
linear combination of the order Q and degree D, while a naive and obvious upper 
bound for N is (Q -I- 1)(D -I- 1). We denote by / the difference between this naive 
upper bound and the actual number of required coefficients N. 

We have no proof of this formula, but it has been found to work 7 for all the 
cases we have considered 

A^== d-Q +q- D -C ^ (g + l)(L> + l) -/. (26) 

This ODE formula is revisited in greater detail in |Appendix B| where all its parameters 
have now been understood. In all cases we have investigated, the parameter q 
appearing in ((26|l is the minimal order of the linear ODE that annihilates S{x). The 

f Of course the minimal order operator right-divides all these Lqu- 

§ The sizes (order and degree) of minimal order vs. optimal ODEs are very well understood in some 
particular cases. For instance, the minimal order ODE (also called "differential resolvent" ) satisfied 
by an algebraic function has coefficients whose degree is cubic in the degree of the function, while 
there exists a linear differential equation of order linear in the degree whose coefficients are only of 
quadratic degree [f4| . To our knowledge, analogous estimates do not exist yet for the (more general) 
case of linear ODEs satisfied by integrals of algebraic functions, such as x'"'- 
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parameter d is the number of singularities (counted with multipUcity) excluding any 
apparent singularities and the "true" singular point x = 0, which is already taken care 
of by the use of the (so-called homogeneity) differential operator x-^. Finally we note 
that the degree of the apparent polynomial of the minimal order linear ODE (as well 
as the other parameters d, q, C and /) can be extracted exactly without obtaining 
the minimal order linear ODE (see (|B.4|) in [Appendix B[ ) . 

As stated above we are dealing with linear differential operators of higher orders 
than the minimal order and whose coefficients are known modulo a prime. To factorize 
such large order linear differential operators, we make use of the method sketched in 
Section m This is done by "following" the series pertinent to a specific local exponent 
at a given singular point. Linear combinations of series with different local exponents 
can be studied as well. Our approach is similar to the one proposed by van Hoeij in ff 5J 
(a "local" factorization deduced from formal series analysis around each singularity 
followed by a "Hermite-Pade approximation" to obtain the "global" factorization). 
The main difference i^|3 that in our case the operators to be factorized are defined 
over a field of positive characteristic 

Actually, the modular nature of our calculation is of great help in this since, with 
the coefficients being known modulo a prime, the coefficients in the linear combination 
of solutions with given local exponents can take only a finite number of integer values, 
so that "guessing" the correct combination can be done by exhaustive search. For 
each series used as a candidate to "break" the linear differential operator under 
consideration we compute three (or more) linear ODEs and from the ODE formula ([26|) 
we infer the minimal order. 

Another point that we address in this paper is the "complexity" of the linear 
differential operators corresponding to the x^"'* seen through the various factors 
occurring in the factorization. One notices that for x^^\ the factors are either of 
order one, or are symmetric powers of the linear differential operator Le- In contrast, 
the linear differential operator for x'^^ contains a factor Z2 of order two which is 
not equivalent to the linear differential operator Le- Recently it has been shown [6] 
that the solution of the linear ODE corresponding to Z2 is a hypergeometric function 
(up to a Hauptmodul pull-back) corresponding to a weight- 1 modular form (see [TT]). 
We think it is important to examine whether the linear ODE for x^^'* contains other 
factors, besides various factors equivalent to symmetric powers of L^, such as the 
factor Z2 occurring for x^^'' j that may have a modular form interpretation. 

Finally we wish to emphasise that the linear differential operator L33 is globally 
nilpotenll§| since it corresponds to a linear ODE that annihilates an integral of an 
algebraic integrand ([3]) (it is "derived from geometry", DFG, see [6] and references 
therein). While this paper is not concerned with global nilpotence as such it must 
be emphasized that the nilpotent condition places very severe restrictions on a linear 

fjNote that the DFactor routine from the DEtools Maple package, corresponding to an 
implementation of these local-to- global ideas I15| . fails to factor Lg of x'^'- The method we display 
in Section |4] actually succeeds in finding this factorization. 

^ Note that in principle one could also resort to algorithms specially dedicated to factoring linear 
differential operators modulo a prime p [16| . However, at present these algorithms are far from being 
efficient enough to handle factorizations modulo primes as big as pr = 2^^ — 19. 

§ The mathematical concept of global nilpotence while quite formal is nevertheless easy to explain. 
Firstly, if p is a fixed prime number, then the differential operator L is said to have nilpotent p- 
curvature iff modulo p, it right-divides the pure power T)^°^'^^^^ of the derivation {T>w = d/dw). 
Secondly, L is called "globally nilpotent" if it has nilpotent p-curvature modulo almost all prime 
numbers p (all primes except a finite number of prime). 
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differential operator, and in particular, provides a proper frameworP^ for the existence 
of basis of series solutions modulo primes for the ODEs (see Theorem 4 in |21|). 
Furthermore, we use two important consequences of the global nilpotence of L33. 
Firstly, globally nilpotent operators are necessarily Fuchsian and permit only rational 
local exponents. Secondly, if L is globally nilpotent so is any factor of L. 

3. Working with non-minimal order linear differential operators 

Once one introduces the particular linear combination (jl4p of x^^^ ^-nd x'-'^^' and 
(modulo the prime pr) of x^^\ it is sufficient to focus on the resulting series and its 
linear ODE (with the operator i29)- From the linear ODE for x*-^-* it is straightforward 
to obtain the (linear) recursion for the series coefficients and using the combination 
p4p calculate as many term^ as required for It is thus not difficult to obtain 

minimal order ODEs requiring fewer terms than pr. We continue however, as in [7], 
to work with non-minimal order ODEs for which fewer series terms are needed than 
for the minimal order ODE. In particular, we make use of the ODE formula (see 
[Appendix B[ ) to infer the order and degree of the minimal order ODE. This formula 
also enables us to control the minimum number of series terms necessary to find 
a Fuchsian linear differential equation (of a given order Q and degree D) which 
annihilates the series. In the sequel, when we say that a linear ODE of order q 
has been obtained, we mean that we have obtained at least three non-minimal order 
ODEs and that the minimal order q has been inferred from the ODE formula. 

From the series (|14|) one can obtain many non-minimal order linear ODEs and 
the resulting ODE formula for L29 reads 

N = 68 Q +291) - 706 = (Q + l)(i:' + l) - /. (27) 

Our understanding of the ODE formula (see [Appendix B| and in particular (|B.4|l ) 
enables us to find Dapp = 1169 as the degree of its apparent polynomial for the 
minimal order operator ^29 without actually producing this minimal order operator. 

In [7] we found that there is a simple rational solution of the linear ODE 
corresponding to L^q (and now also £29) which is the square of x^^\ 

whose corresponding linear differential operator we denote Lf . 

In this paper we use linear ODEs that are not of minimal order to represent 
the minimal order linear ODE annihilating a given series. We also compute the local 
exponents at various singular points of the non-minimal linear ODE and consider them 
as local exponents of the minimal order linear ODE. A remark is in order here. The 
local exponents at w ^ of the linear ODE corresponding to L29 ar(|3 

p = l^ 2^*, 3^ 43, 5^ 6^ 7^, s, 9^ lo, 12. (29) 

In our computation, the non-minimal linear ODE that represents L29 is of order 51. 
One should then really obtain 51 local exponents. It so happens that the 22 "spurious" 

§ The Cauchy-Peano theorem 1201 . which guarantees the existence of series-solutions in the classical 
study of ODEs, does not apply to linear differential equations in positive characteristicl This means 
that in general, a linear differential equation considered modulo a prime number p does not admit a 
basis of power series solutions modulo p, even at an ordinary point. 

t Note that the number of coefficients of all the series used in our calculations does not exceed the 
value of the prime p^ . 

I Throughout the paper the multiplicity of an exponent is denoted by a power: 2, 2, 2, 2 — > 2*. 
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exponents appear in the indicial equation as roots of polynomials in p of degree two 
and higher. These exponents are not rational (indicial polynomials modulo a prime 
of degree higher than one and irreducible) and therefore cannot be local exponents 
for L29, which is globally nilpotent and, hence, has only rational exponents 6 . Had 
the indicial equation of the non-minimal linear ODE given more than 29 rational 
exponents then we would have had to produce other non-minimal linear ODEs and 
obtain the local exponents of L29 as those common to all the non-minimal ODEs. 



4. Factorization of differential operators versus local exponents 

Let us start by giving a brief overview of our factorization procedure. For a linear 
ODE of order q (which may be known exactly or modulo a prime) we compute the 
local exponents at a given singular point w — Ws (such as the origin). We create a 
series Sp{w) starting with the highest integer exponent Up (we seek mainly to utilise 
only those solutions analytic at the origin). This series can be obtained to arbitrary 
length (though shorter than the prime in use) in linear time since we have the linear 
ODE and, hence, the recursion for the series coefficients. We then check to see whether 
the particular solution Sp{'w) is the solution of a linear ODE of order less than q. If 
so, the procedure is repeated for each new factor in turn. If not, we generate the series 
iS'p_i(u') starting at the second highest exponent np_i. The series Sp-i{w) contains, 
via a free parameter, a linear combination of the solution Sp{w). We then let the free 
coefficient of the linear combination vary over the whole (finite) interval given 
by the prime Pr we are using, until a linear ODE of order less than q (if such an ODE 
exists) is found. And then the procedure is repeated. 

For a linear ODE of order q, let Lq denote the corresponding differential operator. 
Consider a singular point w — Wg (for instance Ws ~ 0) and assume the local 
exponents at this point are 

E™, = q, (30) 



mi mo 

Pl , P2 ^ 



where mj is the multiplicity of the exponent pj . In our cases the exponents are either 
integers or rational numbers. Here we utilize only solutions, which are analytic at 
the singular point Wg- So in what follows we consider only integer exponents and we 
denote these as Up. We can then plug the series 

Sp{w) = u;"" + akW^, (31) 

/c>np + l 

into the linear ODE. Demanding Lq {Sp{w)) = will fix all coefficients a^. By 
producing enough terms we can find the linear ODE for the particular series solution 
Sp{w), which is by construction a solution of Lq. The resulting ODE will either have 
order q or order qi < q. In the first case this could mean that Lq is irreducible, or 
Lq does factorize but the factor "responsible" for annihilating the solution Sp{w) is a 
left-most factor. In the second case we have the factorization 

Lq ~ Lq^q^ ■ Lq^. (32) 

To summarise, the series corresponding to the highest local exponent leads to either 
the full ODE or to a "breaking" of the original ODE. 

If the series Sp{w) (corresponding to the highest local exponent) reproduces the 
full linear ODE we turn to the second highest exponent Up-i. In this case, a series 
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starting as w"'''^^ + ■ ■ plugged into the original linear ODE, will yield the expansion 

rip — l 

Sp-i{w) ^ w""-! + flfe w'' + a„^ w"" + Ckw'' (33) 

/c>np_i + l k>np + l 

where all Ofc up to (but not including) are fixed and the Cfc's depend linearly on 
the free coefficient a„p, i.e. 5p_i(w) is a one-parameter solution. The series Sp-i{w) 
is a sum of a series starting as w"''^^ + ■ ■ ■ and the serie^ a„p Sp{w). For generic 
values of the coefficient a„p the series Sp^i{w) will give rise to the full linear ODE. 
But for some values of the coefficients a„ , the series Sp-i{'w) may be the solution 
of a linear ODE of order less than q. This is what leads to the factorization of Lg. 
Intuitively we may hope that such a procedure can work for the following reason. If 
the original operator has many smaller factors this would indicate that there is a basis 
of solutions much simpler than those requiring the full ODE. We don't know this basis 
but by taking a linear combination of two formal "full" solutions (which obviously are 
linear combinations of the basis solutions) it is possible that we can find values of the 
combination coefficients such that the resulting series is a solution of a simpler ODE 
(for these special values some of the basis solutions from the two formal solutions 
cancel each other). 

Similarly, the series solution of Lq that starts at the third highest exponent np_2 
will be a two parameters solution (for simplicity we assume that the exponents differ 
by one) 

Sp^2{w) ^ w'^"-^- +anp_,w"--' +anpW"- + ^ Ckw'^ (34) 

fc>np+l 

where the Ck depend linearly on both the free coefficients a„p_j and a„p. 

To demonstrate how the procedure works in practice we consider the seventh 
order linear ODE for x'--"^^ [THl HI] (denoted iy) . At the singular point w = 0, the 
local exponents are 

l^ 2^, 3, 9. (35) 

Acting with the linear ODE for x'^^ on the series that starts as + ■ ■ (i-e. with 
the highest exponent) 

Sq{w) = + ^ afcu;^ (36) 

fc>10 

fixes all the coefficients. We thus obtain the expansion at w = of x'-'^V^j leading 
to the full linear ODE. Of course, this is not surprising, since the series for x'^'^^ used 
to "generate" the linear differential operator L7 starts as as per ([5]), so the unique 
series Sg{w) must be proportional to x^'^^- The series Sq{'w) cannot be used to "break" 
Lg, since this is the minimal order operator annihilating x^'^'- 

Consider next a series that starts as w"^ +• • •, i.e. with the second highest exponent 

S{w) = +^afcu;^ (37) 

fc>4 

We insert this series into the exact ODE for x'-'^'' and then we solve (term by term) 
the equations arising from iy (S'3(w)) = 0. Doing this we find that the coefficients 04, 

§ Alternatively we can view this procedure as looking at a linear combination of the two formal series 
solutions starting as iii"p-i + • ■ ■ and to"p + • • ■ respectively. 
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as, • • •, flg have to be fixed while the coefficient ag remains undetermined and hence 
enters the series as a free parameter. The remaining coefficients are all given in terms 
of ag. 

S{w) ^ + + 22w^ + 7Aw^ + 417w^ + 1465w® + agw^ (38) 
+ 26839wi" + (36ag - 139067) • + (4a9 + 443325) ■ w^^ + ■ ■ ■ 

The terms in S'(w) in front of the free coefficient ag are the coefficients of the 
series Sg{w). We define S3{w) to be the series obtained from S{w) by setting 
ag = 0. In order to break the operator Lj we look at linear combinations Sa{w) = 
5*3 (w) +aSg{w). For generic values of a the series Sa (w) is annihilated only by the 
full ODE of order seven. However, it is possible that for special values of a the series 
Sa{w) is the solution of a linear ODE of order less than seven. 

We do not know if the "splitting" values of a can be obtained except by a 
"guessing" procedure. The use of modular calculations is very useful in the search 
for the special splitting values. The series Sq{w) and S3{w) can be obtained modulo 
any prime Pr and in the modular calculations a can take its value only in the finite 
range [l,Pr]- If a rational splitting value of a exists it can be found by looking for an 
underlying ODE of order less than 7 annihilating the series S'a(w). In the search we 
use the optimal ODE, which is of order 10 and degree 19 with N = 213. We used the 
prime Pr = 32749 in our search. For each value of a G [1,32749] we calculated the 
series modulo Pr and then looked for an annihilating ODE of order 10 and degree 19. 
For any value of a such an ODE exists and for almost all values N — 213. However, 
for the special values a = 7463 and 7467 we have N = 140 and 206, respectively. The 
decrease in iV is a sure sign that a simpler ODE annihilates Sa{w). In this particular 
case we find that the ODE for a — 7463 is of order four while for a = 7467 the ODE 
is of order six. 

In the case a = 7463 the linear differential operator is Xi • Z2 • A^i = B2-0i-Ni = 
B2 ■ Ti ■ Li, while in the case a = 7467 the linear differential operator is Y3 ■ Z2 ■ Ni. 
These linear differential operators are factors of Lj that were already found in jT3] 
(the indices indicate the orders of the corresponding linear differential operators) 

L7 = Ml ■Y3-Z2-Ni ^ B3-Xi-Z2-Ni 

= B3 • B2 • Oi • iVi = B3 • B2 • Ti • Li (39) 

Remark: We note that the method is not special to formal series. A fortiori, 
it applies to linear combinations of solutions suspected of being parts of a direct sum. 
For instance, removing the series ax*-^-* from the series (modulo a prime) x*^^-' — ^x''^' 
(see p^ ). and, letting a vary in the interval [l,Pr] (recall that a which is a rational 
number appears as an integer modulo a prime), will give (for one value of a) a linear 
ODE of order 29 if the linear differential operator Li for x'^' is in a direct sum in L30, 
the linear differential operator for x^^-* — ^X^^"* • If ^1 h^-d ^ot been part of a direct 
sum the outcome would have been an order 30 linear ODE for all values of a. 

5. Factorization modulo a prime of the linear differential operator L29 

We turn now to the factorization of L29 for which we know that Lf © {Z2 • A^i) is a 
factor. We focus solely on the analytical solutions at w = and we first produce the 
unique series that starts as Si2{w) = w^^ + ■ • ■, where the coefficients in 5*12 (w) are 
given by L2g{Si2{w)) — 0. We found that 6*12 (w) is the solution of an order nine 
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linear ODE (with linear differential operator Lg) with ODE formula: 

N ^ 18 0+ 91) - 73 = {Q + 1){D + 1) - /. (40) 

We know that both (x*-^-*)^ and the solutions of the linear ODE corresponding to 
Z2 ■ Ni occurring in Lj, should be in L29. By explicit checking we found that only Lf 
and Ni are factors of Lg. We can then add the solutions of Z2 ■ Ni to the solutions of 
Lg to produce an 11th order linear ODE (denoted Lu). At the operator level this is 
done by a direct-sum construction Lu ~ Lg® {Z2 -Ni). At the series level used in the 
ODE search programs, it can be done by creating a "generic" solution of Lg (generic 
means a series that gives the full ODE) and then forming a linear combination with a 
generic solution of Z2 ■ Ni to produce a series which is a generic solution of Lu. The 
resulting linear differential operator Lu has the ODE formula: 

24 + Hi? -Ill = {Q + l){D + l) - f. (41) 

We have thus shown the following factorization of L29 

^29 = Lis ■ Lu. (42) 

Before proceeding we wish to clarify the meaning of the ODE corresponding to 
the left factors. Having obtained the linear differential operators L29 and Lu, a right 
division should give the linear differential operator Lig. One should bear in mind that 
the order of these operators is large and our representation of them are not of minimal 
order. In the computation, the linear differential operators representing L29, Lig and 
Lii are of orders 51, 32 and 17, respectively. Our representation of the factorization 
(|42|) reads in fact 

022- L29 = (Oi4-Li8)-(06-ill) (43) 

where the equality stands for "both sides acting on S{w) give zero" . Since the series 
solution S{w) demands an order 29 linear ODE, and the order of the extra operator 
Og is arbitrary, there are intertwinners leading to 

O22 • ^29 = Oi4 • da ■ (Lis • ill) (44) 

With the relation Lig ■ Oq — Oq ■ Lig, the linear differential operators Lig and Lig 
are equivalent and have the same factorization structure. 

Next we take the series S{w) = x^^^ — X^'^V^ + x'^Vl-^O and plug it into the 
linear ODE for Lu to produce a new series whose linear ODE corresponds to the 
linear differential operator Lig. This linear ODE has the formula: 

N ^ UQ +18D +873 = iQ + l){D + l) - /. (45) 

To proceed further with the factorization, we compute the local exponents a,t w = 
for the linear ODE corresponding to Lig. These are: 

p = l^ 2, 3^ 4, 5^ 6^ 7^ 8, 9^ lO. (46) 

For the linear differential operator Lig we look at the solution that starts as Sio(w) = 
w^^ + • • • Unfortunately, this gives linear ODEs with the same ODE formula as Lig, 
that is the series reproduces the complete linear ODE represented by Lig. This means 
that the factor responsible for the annihilation of Sio{w) occurs at the left of Lig. 

The second highest exponent is p = 9. When Lig is applied to a generic series 
Sg{w) = + ■ ■ ■ we obtain a series that depends on the coefficient in front of w^'^. 
This one-parameter series starts as (modulo the prime pr) 

+ aio w 1" + (15419 aio + 10040) • u;" + • • • (47) 
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The series, collected in oio, enables us to reconstruct the full linear differential operator 
Lis, but it is possible that the above combination may give a linear operator of smaller 
order for particular values of aio, which have to be found by experimentation. It is here 
that the modular calculations allow us to find a definitive answer. Modulo the prime 
Pr, the coefficient aio spans a finite set of integer values [l,Pr]- The determination of 
the coefficient aio is thus feasible in a finite computational time by exhaustive search. 

With the value aio = 12999 we found that the series Sg gives a linear ODE of 
smaller order than Lis with ODE formula 

N= 36Q +13£) +715 = (Q + + 1) - /. (48) 

The local exponents at w = for this linear ODE of order 13 (that we denote Lis) 
are: 

p = 1^ 2, 3, 4, 5^, 6^, 72, 8, 9. (49) 
The highest exponent is indeed nine, of which the associated solution gave us Z/13. 



Let us be more explicit on the meaning of the combination + 12999 
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that has given rise to the 13th order linear ODE. By "following" the series H 

we obtained the full linear ODE. The factor that annihilates this series is thus to 
the left in the factorization of iis. We find that this series comes with a log(w)'* 
behaviour (see below). By the combination w'^ + 12999 w^'^ + • • •, we are looking for 
the particular value of aio that removes this logarithmic solution from the linear ODE 
corresponding to Lig. Thus the linear ODE for L13 no longer has a solution behaving 
as log(u')'*. 

Completing from L13 to Lig we obtain a fifth order linear ODE (called L5) with 
ODE formula: 

N= 8Q +5£) +912 = (Q + 1)(£) + 1) - /. (50) 
We thus have the factorization: 

L29 = ^5 • -^13 • ^11- (51) 

The fifth order linear differential operator L5 is the one whose existence we conjectured 
previously and which we believe should annihilate a homogeneous polynomial of the 
complete elliptic integrals E and K of (homogeneous) degree four. 

The local exponents at the origin of the linear ODE corresponding to L5 are 

w= 0, p = 1, 3, 6, 9, 10. (52) 

Plugging a generic series ^ c„ into the linear ODE fixes all the coefficients 
(including ci, C3, cg, cg) with the exception of the coefficient ciq. The "survival" 
of a single coefficient is a particular feature of an irreducible factor with one non- 
logarithmic solution. The exponents at the other singularities (apart from w = 00) 
are: 

w = 1/4, p= -292, -28, -23, 0, 

w = -1/4, p= -352, -33, -31, 0. 

This suggests that one should plug the following ansatz into the linear ODE: 

4 

' P4-i,i • K^-' (53) 



(l-4«;)29(l + 4i(;)35 ^ 



i=0 

where K and E denote the complete elliptic integrals 

K = 2Fi([l/2,l/2],[l],16u;2), E = 2Fi([l/2, -1/2], [1], W). 
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Collecting terms of the form K^~''^ E'' we determine the polynomials P^-i.i whose 
degrees are increased steadily until we obtain a solutior^. With degree around 200 
the following solution was found 

^(1 - IQui^ f F4,o ■K'^ + (1 - 16^2)^ P3,i • E 



(l-4w)29(l + 4u;)35 

+ {1 - \Qw^) P2,2 ■ E^ +Pi,3-KE^ +PQ^i-E* 

where P^-i^i are polynomials in w with coefficients known modulo the prime pr and of 
degree respectively 200, 202, 204, 204 and 204. The expressions for the polynomials 
Pi-i^i can be found in (22j . As conjectured the linear differential operator L5 is thus 
equivalent to the symmetric fourth power of Le- 

5.1. The linear differential operator Ln has six factors 

As shown above the linear differential operator L29 factorizes into three factors of 
order 11, 13 and 5. We have just shown that the fifth order linear differential operator 
is irreducible. Next we consider the linear differential operator Ln. 

We know that the fourth order linear differential operator L\ © (Z2 • iVi) is a 
right-most factor of Ln, so we obtain 

Lii - Nt{L{®{Z2-Nx)). (54) 
The ODE formula for the seventh order linear differential operator N-; reads: 

N= 15 0+7D+89 = (0 + l)(i:' + l) -/. (55) 
At It; = 0, the local exponents (for iVy) are 

p = 2^, 3, 42, 5, 12. (56) 

Plugging the series +X]fc>6 ^k'W^ into the linear ODE for TVV fixes all the coefficients 
except ai2 corresponding to a solution with the highest local exponent. Letting the 
combination coefficient 012 vary in the finite range [l,Pr-], we found for the value 
ai2 — 22292 a linear ODE of order less than seven, with ODE formula 

N= HQ + bD +1% = {Q + 1){D + 1) - f, (57) 

and with exponents at the origin 

p = 22, 3, 4, 5. (58) 

For this linear ODE we consider the one-parameter series that starts with w'^ and 
which contains the coefficient 05 as a "free" parameter. By letting the coefficient 
05 vary in the finite range [l,Pr], we found that for 05 — 29103, the linear ODE of 
order five breaks into two linear differential operators of order three and two, O3 • O2, 
respectively. 

In conclusion we have decomposed the differential operator Ln of order 11 into 
six irreducible factors 

ill = O2 • O3 • O2 • [LI ® (Z2 • iVi)), (59) 

where the indices denote the order of the corresponding linear differential operators. 



^ Once the solution has been obtained a check to any order can be carried out. Typically a good 
check amounts to plugging polynomials of degree 300 into (I53I I. 
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5.2. The linear differential operator Ln in exact arithmetic 

B. Nickel has obtaine4i some linear differential operators that right-divide and, 
especially, the linear differential operators (that we call) U2 ■ Ni and F3 ■ F2 - if. 
Checking these operators against the factorization ([SQ]) . we found that Ln has the 
direct-sum decomposition 

ill = {Z2-N,) ®V2(B (Fs-F^-Ll), (60) 

where V2 is equivalent to U2 (or to O2) and the product F3 ■ F2 is equivalent to the 
product O3 • O2 in ([551) . 

Furthermore, using some tricks in the modular calculations supplemented with 
constraints on the apparent polynomials ([13], Appendix A in [7]), Nickel succeeded 
in finding the considered linear differential operators exactly. The linear differential 
operators V2 , F2 and -F3 are giverQ in [Appendix C[ 

The procedure for obtaining a rational number from its image modulo a prime 
is known as "rational reconstruction^^ and has many applications (for details see 
e.g. [Ml [2SI IMI)- Consider a rational number n/d which has the residue u modulo 
the prime m. Given u and m, a rational reconstruction algorithm tries to recover the 
rational n/d under some conditions on the magnitude of the unknowns n and d. The 
simple version of this condition is 

2N^ < m, N = max{\n\,d). (61) 

The algorithm will then output the rational n/d satisfying the above condition, but 
this rational number may not be the actual one for the problem. If the residue is 
known for several primes mi, it is the Chinese remainder u which is considered and 
m is the product of the primes rui. 

In any case, the knowledge of the order of magnitude of n and d is important. For 
instance, the exact linear differential operator F2 can be recovered using the results for 
three prime^ In our modular calculations the residues are coefficients of polynomials 
occurring in Fuchsian linear ODEs. Besides the order of magnitude, which can be 
guessed, the linear ODE, once reconstructed, should satisfy certain properties. The 
indicial equation should give the correct local exponents, which have been obtained 
either from a linear ODE modulo a prime or from a diff-Pade analysis [7]. For the 
apparent singularities the conditions on the apparent polynomial should be verified. 

Consider the linear differential operator F2 (see [Appendix G\ 

F2 - V2{w) Papp{w) ■ Dl + Vi{w) Pi{w) ■ D^ + Po{w). (62) 

The singularities are known and are roots of the polynomial 7^2 (w) (the polynomial 
Vi{w) contains a subset of the singularities). We note that Papp{w) (the apparent 
polynomial) and Pi(w) can be reconstructed with two primes while Pq^w) demands 
three primc^. 

However, when Papp{w) and Pi{w) have been found, the matching of the known 
local exponents will fix some coefficients in Pq^w), but, more importantly, we get 
an idea about the order of magnitude of the common denominator in the various 
rational coefficients in Pq{w). This magnitude was found to be 2^^ and, with this 

tt Private communication. 

t We are grateful to B. Nickel for these results. 

J We have actually obtained the linear differential operator L29 = Lig- Ln for four primes 2^^^ — 19, 

- 49, - 51 and 2i5 - 55. 
§ We have in each case one or more additional results modulo a prime for checking our calculations. 
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scaling, the still unreconstructed coefficients in Pq(w) can be reconstructed using two 
primes and then checked against the conditions for the apparent singularities. Our 
iterative procedure for reconstructing the exact polynomials thus amounts to first 
reconstructing the polynomials Pappiw) and Pi (w) (with two primes) in order to obtain 
Po{w) with two primes instead of three. More details can be found in [Appendix D[ 
which deals with a rational reconstruction experiment on the apparent polynomial of 
F3. One should note that, to rationally reconstruct a linear ODE, it must be obtained 
with as many primes as necessary. When the results are time consuming or hard to 
obtain, the knowledge of the underlying problem may help one to guess the scaling 
factor which forces the condition ((6T|l . 

Note that when the number of primes is not sufficient some reconstructed 
coefficients will, obviously, be in error. A strong check can finally be done on these 
linear differential operators that should convince one of their correctness. The linear 
differential operators F2 and ^3, that we are looking for in exact arithmetic, are factors 
of the linear differential operator L33. They are necessarily globally nilpotent [Bj since 
I/33 is. The linear differential operator L33 is globally nilpotent since it corresponds to 
a linear ODE that annihilates an integral of an algebraic integrand ([3|) (it is "derived 
from geometry", DFG, see [6] and references therein). We have calculated the p- 
curvatures of these reconstructed linear differential operators F2 and F3, and found 
that they are, indeed, globally nilpotent. 

The global nilpotence of a linear differential operator is a strong and very special 
property that is rigid enough to make us totally confident that the polynomials 
occurring in the linear differential operators F2 and F3 have been reconstructed 
correctly. 

As for the solutions of the factors occurring in Ln, the simple V2 is equivalent to 
the linear differential operator Le and its linear ODE annihilates 

(63) 



1 — 4:W V 1 — 16w 
From ((63)) it is straightforward to see that V2 is actually equivalent to the second order 
operator for x^^^^ (denoted Nq in [H]). Note that V2, equivalent to A'o, does not mean 
that x^^"* itself is a solution of L29, but rather some linear combination of x*^^' and 
its first derivative is. Indeed ([551) can be expressed as: 

(l + 4^)(l + 8^.^) dx(^ -4.(1 + 4..).X(^). 
2w dw 

The remarkably simple result ((63)) begs the question about the very nature of V2. Is 
the occurrence of x'^'' (and its first derivative) in x'"^'' a mere coincidence or does it 
suggest a more general structure. Does an operator equivalent to x^'*'' appear in x'^' 
(or x'-'^-' in X*-^-*)? If we do not have this very strong result, is it nevertheless the case, 
that some of the individual factors occurring in say the factorizations of L7 appears 
in x^^''? Expressed more generally is it the case that operators equivalent to factors 
from x'™-* appear in the factorization of x*^"^ (with m < n and n and m of different 
parity)? Similar questions can be asked with regard to the occurrence of the rational 
solution of Lf that we have written as (x''^'') in (|28p . We have already conjectured in 
(|23p that for the x'-"-''s we have direct-sum structures corresponding to linear ODE's 
of smaller order for selected linear combinations involving x^"^^'^'' . Could it be that 
polynomial (i.e. non-linear) combinations of x'''"^'s {m < n) can be used to further 
deplete x*"""* or at least appear as factors? 
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As emphasised in [B], we have a strong belief (but no proof) that aU the Unear 
differential operators occurring as factors of the hnear differential operators for the 
^("^s, have to be related to the theory of elliptic curves (complete elliptic integrals E 
and if, algebraic modular functions expressed as algebraic hypergeometric functions, 
modular forms of weight-1, etc.). The calculation of the p-curvature of F2 and F3 
excludes linear differential operators associated with algebraic functions. The simple 
occurrence in of E and K is in contrast to the linear differential operators F2 
and i<3, which do not seem to be equivalent to a symmetric power of Le- We must 
explore whether (similarly to what we found [5] for Z2 and the linear differential 
operator occurring in the analvsi jjfl of <i>^'') the linear differential operators F2 (resp. 
F3) correspond to modular forms of weight-1 (or higher weights), or 2F1 (resp. 3F2) 
hypergeometric functions with a Hauptmodul pull-back (up to multiplication by some 
n-th root of a rational function). 

To see if a solution of the second order operator F2 is a 2F1 hypergeometric 
function with a Hauptmodul pull-back (up to multiplication by the n-th root of some 
rational function) would require one to find not only the Hauptmodul pull-back, but 
also a change of variable (covering) mapping the large set of singularities in F2 onto 
three singularities (0, 1, 00), and find, besides, the rational function occurring in the 
multiplicative factor in front of the hypergeometric function 2F1 (which looks like the 
Hauptmodul [5]). 

The occurrence of an involved apparent polynomial is a quite severe obstruction 
for performing these educated guesses. It is always possible to get rid of the apparent 
polynomial of a linear differential operator by introducing a higher order, but still 
Fuchsian, linear differential operator with no apparent singularities (see [Appendix C] ). 
This however is not helpful. What we really need is to find an equivalent linear 
differential operator with a smaller apparent polynomial or, hopefully, no apparent 
polynomial. 

This is how we achieved [6] such a calculation for the linear differential 
operator Z2. We were able to find a second order operator (occurring as a factor 

(3^ 

in ), which is simpler than Z2 because its apparent polynomial is just 1 — 2w 

Z2-Mi = Mi-Z2, (64) 

where Mi and Mi are two first order linear differential intertwiners. Up to the change 
of variable (covering) 

^ " Ti VFT — r^' 

(1 — w) (1 — iw) 

wrapping the seven singularities of Z2 onto the three singularities of the 
hypergeometric function, we were able to find a modular form of weight-1 solution 
of the equivalent second order operator Z2. This was a consequence of its very simple 
apparent polynomial. At the present moment, we have not been able to replace F2 by 
an equivalent second order operator with a simpler apparent polynomial. The situation 
is even more involved for F3 (see [Appendix C| ). The modular form interpretation of 
F2 and F^ remains to be done and is clearly a worthy challenge. 

tfThe functions '3?^' arc simplified Ising type integrals [9] obtained from the x^"' integral 
representation JSj by setting the Fcrmionic factor G'"' = 1. 
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5.3. The linear differential operator L13 

We continue our procedure of factorization for the 13th order operator L13 occurring 
as a factor in L2g. Recalhng the local exponents at w = 

p = 1^, 2, 3, 4, 5^ 6^, 7^ 8, 9 (66) 

we know that the series corresponding to the highest exponent enables one to 
reconstruct the full linear ODE. The next series to consider is thus the one-parameter 
series + agw^ + ■ ■ ■. For every value of the linear combination coefficient og in the 
interval [IjPr-]; we found that the resulting linear ODE is of order thirteen. 

Both series {w^ + ■ ■ ■ and + agw^ + • • •) are annihilated by a left-most factor 
in L13. To proceed with the factorization, we would have to consider "deeper" 
combinations of series solutions (see Section 2]). For instance a.t w — 0, we could 
use the two-parameter (ag, og) solution + asw^ + agw^ + • • then the three- 
parameter (07, as, ag) solution + arw"^ -I- as + agw^ + • • and finally the 
four-parameter solution + aew^ + arw"^ -I- as + agw^ + ■ ■ ■. However, if to 
is the computational time for a single ODE search, then to check the factorization 
using a solution with k free parameters requires a computational time of Pr ■ This 
requires a very long time for the prime Pr = 32749 taking into account the sizes of the 
linear ODEs we are dealing with here, and hence we have not pursued this approach 
beyond the one-parameter case. We could also use, from the outset, the most general 
five-parameter solution +aQW^ +a'jw'^ -I- as -I- ag + • • •> (only as = and 

as = 1 need to be considered) that should give all the factorizations (if any) of L13. 
However, a check of the factorization using this five-parameter series solution clearly 
suffers from the prohibitive time requirements mentioned above and it is beyond our 
current computational resources. 

In Section|4]we described our method of factorization modulo a prime by focusing 
on the singularity at the origin. This singular point has no special properties 
which makes it better suited than other singular points for our factorization scheme. 
However, to have a clear working scheme, the singular point one chooses to focus on 
must be sufficiently singular (by which we mean, in this case, have many confluent 
logarithms) to allow one to extend from the local scheme to the global scheme. So, we 
looked at what happens if we use expansions about other singular points. 

Considering the ODE corresponding to L13 translated tc0 w = 00 one can 
follow the series of the highest exponent which is —30. This series also demands 
the full ODE. The one-parameter series corresponding to the second highest exponent 
x~^^ 4- a-30 x~^^ + • • •, (with x — l/w) also gives rise to the full ODE (i.e. the order 
remains 13) for all values of a_3o S [l,Pr]- Similar calculations were performed for the 
linear ODE translated tcl3 ^ — 1/4- Neither the series -(- • • •, nor the one-parameter 
series cc^ + a^^x^ + ■ • (with x = w — 1/4) gives rise to a linear ODE of order less 
than thirteen for any value of 05 g [l,Pr]- The series solutions in front of the higher 
logarithmic solutions around some other singular points give rise to the full linear 
ODE of order thirteen. 

Instead of considering the series solution as +a%'w^ -\- 07 -I- as ui^ -I- ag -\- 
• • •, with its prohibitive computational time requirements, we decided to try another 
procedure that may give us an idea about the number and order of factors occurring 
in ii3 (if reducible). 

1 The local exponents at ui = 00 are -40, -39^, -38, -36^, -35, -34^, -33, -32, -31, -30. 
tt The local exponents at ui = 1/4 are -9, -8, -7, -6, -5, -4, -3, 0, 1, 2, 3, 4, 5. 
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We start by examining how the various formal solutions of L13 appear. Consider 
(near w = 0) a, general solution with log's such as 

S{w) = Sn{w) logH" + Sn-i{w) log(«;)"-i + • • • + So{w), (67) 

where the exponent n of the log is generically taken to the maximum allowed value of 
12, i.e. the order of L13 minus one, and where the Sn{w) are power series expansions 
Oq"^ + a^"^ w + ■ ■ ■. Plugging the solution S{'w) into the linear ODE corresponding 
to ii3 and solving I/13 {S{w)) = term by term, we found that the highest allowed 
exponent is n = 3. 

We therefore fix the solution S{w) as 

S{w) = Ssiw) log{wf + S2{w) log{wf + Si{w) log{w) +So{w), (68) 

and act on it by the linear ODE corresponding to L13 and solve term by term (we have 

to solve only to since this is the highest local exponent for L13 around w = 0). 
The coefficients (up to w^) in the Sk{w) must be fixed. Among the 40 coefficients 13 
will remain as free parameters (equal to the order of the linear ODE). Attached to 
any of these free coefficients is an independent solution of the linear ODE. 

To clarify the scheme of these solutions, from which we shall infer the number of 
factors, we solve L13 {S{w)) = 0. This leads to the equation 

E(EC'(")log(«;)")«;'= =0, (69) 

fe>0 n=0 

which we solve for each k and n by using the following recipe: The coefficient Cj^^ 
of the term w'' log(tt;)" will in general be a linear combination of coefficients a^™^ 
from Smiw) with j < k and n < m. When the coefficient c|,"'' contains coefficients 
aj™'' with m = n only, we solve for the coefficient aj"'' of highest index j. When 
the coefficient C^"'' contains coefficients Oj™'' with to < n, we solve for the coefficient 
aj™'' of lowest index m and highest index j. This is because we know, from all our 

computations on Ising type ODEs, that when a solution such as S{w) log(w)("'' H 

occurs, S{w) log(w7)'^"^-') + • • (with the same S{w)) is also a solution. 

We introduce the notation [w^] to mean that the series begins as {const. -{ ). 

The results of the computation are the following. Four solutions can be written as 

[w"^] log{wf + [w^] log(w)^ + [w] log{w) + [w], 
[w"^] log(w)2 + [w^] log(u;) + [w], 

[w"^] log(w) + [w], and [w^]. (70) 

Four other solutions can be written as 

[w^] log(w)^ + [w^] log(w7)^ + [w] \og{w) + [w], 
[w^] log(w)2 + [w^] log(w) + [w], 



log{w) + [w], and [w^]. (71) 

There are also two sets of two solutions each that read 

[w^] log(w) + [w], and [w% (72) 

and 

[w^] log(w) + [w], and [w% (73) 
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Finally there is a non- logarithmic solution starting as + • • •. 

In view of this scheme, one may conclude, in analogy with all the Ising calculations 
we have performed and where hypergeometric functions occur, that the factors 
occurring in L13 are of order 4, 4, 2, 2 and 1. At the point w = 00, one obtains 
the same structure of solutions leading to the same scheme, that is factors of order 
4, 4, 2, 2 and 1. However, at the singular point w = 1/4 the structure changes slightly. 
The solutions, grouped as done above for the point w = 0, lead to a scheme of six 
factors with orders 4, 2, 2 2, 2 and 1. 

To reconcile the threqjj schemes (around w = 0, w = l/i and w = 00), the linear 
differential operator L13 may have either four factors of orders 4, 4, 4 and 1 or five 
factors of orders 4, 4, 2, 2 and 1. 

Around the three singular points, the schemes allow for an order one differential 
operator whose corresponding series starts a.s + ■ ■ ■. It happens that this first order 
differential operator (call it Li) occurs as a right-most factor of L13 = L12 • Li. 

We have found that the solution of the linear ODE corresponding to Li is a simple 
polynomial of degree 34, which reads (modulo the prime PrJiM 

P{w) = u)^ + 30849 + 4080 + 11244 + 26721 w;^ 
+ 29301 + 23070 + 30185 u;^^ 26217 + 10853 
+ 25659 w^^ + 4536 w^^ + 31400 + 22061 w^^ + 31481 w^^ 
+ 3767 w'^° + 6508 w^^ + 10160 w^^ + 31426 w^^ + 29441 w^^ 
+ 17755 + 6024 w^^ + 31840 + 10393 w^* + 20669 w^^ 
+ 4477 + 29192 w^^ + 20075 w^^ + 2957 w^^ + 2003 w^^. (74) 

Although we have obtained such polynomials for the four primes, previously 
mentioned, our attempted rational reconstruction [24l [25l [2^ of the exact P{w) was 
not successful. There is no further information to guide our quest for the exact P{w). 
There are only two indicial exponents (5 and —34) corresponding to the two points 
w — and w = 00, respectively. The linear differential operator ii is a first order 
linear differential operator of the form: 

~ d dln(P(u;)) 

Li ^ — + . (75) 

aw aw 

It is thus automatically globally nilpotent. Global nilpotence is a very severe constraint 
to fulfill for higher order linear differential operators, but for first order operators like 
(|75p it provides no additional constraints on P(w). 



6. Comments and speculations 

In view of the previous results, we give some concluding remarks. The linear 
differential operator L29, corresponding to x^^^ — x^'^V^ + x'^^Vl-^O, can be written 
as 

L29 = ^5 ■ ■ ill, (76) 

with 

Lii = (Z2 • TVi) ® F2 © (^^3 • F2 ■ LI) , (77) 
L13 = L12 ■ Li. (78) 

t There are not enough logarithmic solutions at the other singular points. 

J With the solution P{w), we have the coefficients for the deepest combination series w^' + ag + 
07 ui'' + asw^ + agw^ + ■ • •. 
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The linear differential operator L5 is equivalent to the symmetric fourth power of 
Le- This linear differential operator is the factor of maximum order, assuming that 
the factorization scheme of L12 is correct. The scheme of factorization ([7S)) then 
generalizes what we have obtained for x^^-* and x'^''- 

Our conjecture on the structure of the x^"\ namely, (|23|) and ((24| would give, for 
the six-particle contribution x^^-* , the following scheme 

X^') = a4-x<^) +a2-x'2) (79) 

with a solution of a linear ODE of order q whose corresponding linear differential 
operator factorizes as Lq = Lq ■ Lg^g, and where the left-factor Lg is equivalent to 
the symmetric fifth power oi Le- 

As far as the singularities are concerned, the 11th order linear differential operator 
Lii has only the singularities of the linear ODE corresponding to Lj (the operator 
for x^^^) and the "unknown'!^ w — 1/2. This singularity w = 1/2 occurs in the 
third order linear differential operator F3. The second order differential operator F2 
is responsible for the p = —5/4, p = —7/4 singular behavior around the ferromagnetic 
point w = 1/4 (see Table 4 in |7j). 

The singularities of the linear ODE corresponding to the block L12 are (besides 
w = 0, ±1/4): 

(1 - w){l + 2w){l + 3w + 4u;2) 

(1 + w){l - 3w + w^)(l + 2w- A'u?){l - w - 3w^ + 4w^) 

(1 + 8w + 20w^ + \bw^ + 4w'*)(l - 7w + bw^ - Aw^) 

(l + 4w + 8u;2). (80) 

The singularities in the first line are also singularities of the linear ODE for x'^'- Note 
that at w — —1, w = 1 and w = —1/2, for instance, the linear differential equation 
corresponding to L12 has logarithmic solutions. Therefore, at least one of the factors 
(if the scheme is correct) in the block L12 is not equivalent to a symmetric power of 
Le- If we consider the possibility that the linear differential operator of order twelve 
L12 is irreducible, this would mean that we are faced with a highly restricted object, 
which is globally nilpotent and not eauivalenl Fpi] to the symmetric eleventh power of 
Le. 

Let us close with the following question arising from some of the modular 
calculations and rational reconstructions presented in this paper. Is it in general 
easier to generate series for many primes, use these to reconstruct the exact series 
and hence obtain the exact linear ODE, or is it easier to obtain the linear ODE for 
a smaller number of primes and then carry out the rational reconstruction on the 
coefficients of the linear ODE? To disabuse the reader of the obvious first impression 
that the second method must be easier, we would like to point out that when we opt 
for a non- minimal order linear ODE, we gain by way of a reduction in the number of 
terms necessary to find the linear ODE, but this comes at a cost of an increase in the 
size of the coefficients in the linear ODE (see the last paragraph of Appendix A in [7] 
for an estimate for x*''^^)- Even if we are dealing with the minimal order linear ODE, 
the coefficients in the right-factors have less digits than the coefficients occurring in 
the left-factors. For instance, considering the first factorization in (I39p . the maximum 
number of digits is 6 in Z2, it increases to 27 for ¥3 and to 33 for Mi. 

§ Unknown with respect to the "t"^' integrals [S] and to our Landau singularity analysis [?]■ 

tf But might, for instance, be equivalent to a symmetric power of a smaller order globally nilpotent 

operator related to modular forms. 



High order Fuchsian ODE for x*-^-* 



22 



7. Conclusion 

Using tiie Fuciisian linear ODE of x'^' (obtained modulo a single prime pr), we 
have been able to go quite a long way towards understanding the factorization of 
its corresponding (minimal order) linear differential operator L33. In particular we 
have found several quite remarkable results. 

The direct-sum structure of L33 generalizes what we have found for the linear 
differential operators of x^^"* and x^^^- the linear differential operators of x'^' '^e 
found not only the occurrence of a term proportional to x'^'' j but also the occurrence 
of a term proportional to x''^''- We conjecture that this structure occurs for all x*'"\ 
i.e. we expect to see terms in x'"'' proportional to x^""^*^-*- 

The linear differential operator L29 annihiltating the "depleted" linear 
combination (jl4p for x^^'' follows the same structure seen for x*''^'' a-nd x'^''- The 
left-most factor of L29 is equivalent to the symmetric fourth power of the second order 
operator corresponding to complete elliptic integrals of the first (or second) kind. 

Some right-factors of L29 are given in exact arithmetic. In particular one notes 
the occurrence of a very simple second order operator V2 and of the remarkable factor 
Z2 that occurred in x^'^''- Using the series of x*'^^ obtained for pr as well as three 
additional primes, we have obtained the linear ODE modulo each prime and checked 
that the mentioned right factors are indeed exact. We have used the results for the 
four primes to see whether a rational reconstruction of right factors is feasible. 

Two of the right- factors, F2 and F3 (of order two and three, respectively), are 
highly restricted globally nilpotent linear differential operators, but, unfortunately, we 
have not been able to find exact solutions as we did for the linear differential operator 
Z2 occurring in the factorization of the linear differential operator for x*^'^''- Providing 
a better understanding of these operators, say in terms of modular forms, is clearly a 
natural (but actually quite difficult) challenge. 

The incomplete part of our analysis is concerned with the 13th order linear 
differential operator L13. For this operator we did find a right-factor of order one 
which, quite remarkably, has a polynomial solution. The factorization of the remaining 
12th order linear differential L12 is beyond our current computional ressources. By 
producing all the twelve formal solutions of L12, a factorization scheme appears 
where the differential operator L12 (if reducible) could factorize into three fourth 
order operators, with a possible scenario that one of the fourth order operators could 
factor into two second order operators. Clearly some work remains to be done to 
better understand L12, and hopefully find the exact fourth order operators in its 
factorization. We thus hope to gain a better understanding of their very nature (are 
they symmetric powers of Le ov perhaps linear ODEs associated to modular forms, 
namely hypergeometric functions with a Hauptmodul pull-back) . 
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Appendix A. Solution of the order four linear ODE M2 occurring in x'^^^ 

Defining x ^ and 

K = 2F1 ([1/2, 1/2], [1], 16.T) , E ^2Fi ([1/2, -1/2], [1], 16a;) (A.l) 
the solution (analytical at a; = 0) of the linear ODE corresponding to M2 is 
1 



a;4 • (1 - 16a;)4(l - Ax){7 + 80x) 



(l - 16x) ■ P3,o ■ - 3-P2,i-K'E 



with: 

P3,o = 819200 - 1050624 a;"* + 494976 a;^- 39128x^-90 a; + 63, 

P2.1 = 26214400 + 1458176 - 4698112 a;"* + 678464 a;3- 26120 a;2- 818 x + 63, 

Pi,2 = 7208960 x^ + 1169408 x"* - 300288 + 8728x2 + 538 X- 63, 

Po 3 = 363520 x'' + 53696 x^ - 1144x2 -86 x + 21. 

Appendix B. The ODE formula 

The linear differential equations annihilating a series S{x) we are interested in are 
Fuchsian. This means that all singular points of the linear ODE, and in particular 
X = and x = 00, are regular. A form of the linear differential operator that 
automatically satisfies this constraint is: 

Q D / d \* 

^QD = '^0 V dx ) ' "^30^^' "-Q^ ^ (^--^^ 

i=0 j=0 ^ ^ 

The condition ago 7^ (resp. 7^ 0) is required to make x = (resp. x = 00) 

a regular singular point. 

Note that it is the use of the (homogeneitj{^ operator x ^ (rather than just 
d/dx), which leads to the above conditions guaranteeing the regularity of the singular 
points X = and x = 00 and to the equality of the degrees of the polynomials in front 
of the derivatives. For the operator d/dx, a simple rearrangement of terms shows that 
(|B.1[) can be re- written as 

i=0 j=0 ^ ' 

where the coefficients hij are linear combinations of the . This is the form of the 
Fuchsian linear ODE we have used in many previous papers (e.g. [8| [9l ITO l Fn ] [27 t l28] ) . 
The Fuchsian character of the ODE is reflected in the decreasing degrees of the 
successive polynomials in front of the derivatives. 

To find the linear ODE annihilating S'(x), the coefficients a.y in (jB.ip have to be 
determined. This can be done by demanding Lqd (S'(x)) = 0, resulting in a set of 
linear equations for the unknown coefficients a^. In [7] this set of linear equations 

^ Also called Euler's operator. Recall that [x ;t^)" ■ = k" x'^ . 
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was put into a well defined order and if the corresponding Nqjj x Nqjj determinant 
(with Nq£) — {Q + 1){D + 1)) vanishes, a non-trivial solution exists. The zero- 
determinant condition was checked by creating an upper triangular matrix U using 
standard Gaussian elimination and a solution exists if we find U {N, N) = for some 
A^. The N for which U{N,N) = is thus the minimum number of coefficients 
needed to find the linear ODE for given Q and D. The deviation between the actual 
number of coefficients needed N, and the generic (maximum) {Q + 1){D + 1) was 
called A in [7] . 

To fully understand the deviation A = (Q + 1){D + 1) — N , we may alternatively 
compute the nuUspace of the matrix U. The dimension of the nullspace, if a solution 
exists, is related to A. In others words, solving Lqu {S{x)) — term by term will 
fix all the coefficients but leaves / coefficients unfixed among the Nqd ones. These 
are all independent ODE solutions for given Q and D. 

In [7] we reported a remarkable formula arising from empirical observation 

N ^ d-Q +q- D -C ^ {Q + 1){D + 1) -/, (B.3) 

where we have replaced the parameter A used in [7j by the parameter / that we 
now can understand as the number of independent solutions for given Q and D (this 
understanding will be useful later). The ODE formula (|B.3[) should be understood 
as follows: For a long series S{x) we use three (or more) sets of Q and D and solve 
Lqd i^ix)) = (by nullspace or term by term). From this we obtain the value of 
the parameter / (if / > 0, otherwise we increase Q and/or D) for each pair {Q,D). 
These values {Q, D, /) are then used to determine d, q and C in (|B.3p . In all cases 
we have investigated, the parameter q is the order of the minimal order linear ODE 
that annihilates S{x). The parameter d is the number of singularities (counted with 
multiplicity) excluding any apparent singularities and the "true" singular point x = 
which is already taken care of by the use of the differential operator x-^. 

We revisit in Table IBll some ODE formulae from Table 1 of 7J . We give the 
value of the parameter / corresponding to the same Q and D considered in jTj. The 
first observation is that, generally, both ODE formulae (in Table [BT] and in T) agree. 
When they do not, the difference is in the parameter C. But we remark that C — f 
always equals C — A, which is easily understood from the equality in (jB.Sp . The 
second observation is that, for the linear ODE which have the constant as solution 
(i.e. x''^' a-iid Gx*-^' — Sx^'^-'), the parameter q appears as the actual one. 

Table Bl. ODE formula for x'"', n = 1, 2, ■ ■ • , 5 and for the combinatfons 
g£(n+2) _„^(n)^ n= 1, 2, 3. The last column gives the value of the parameter 
/ corresponding to the same Q and D considered in [7]. 



Series 
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^qD- 
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/ 




IQ 


+ ID 


+ 1 
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1 


1 
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+- 2D 


+ 1 
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1 


1 




12 Q 


+- 7D 


- 37 


11 


17 


2 




7Q + 


10 D 


- 36 


15 


9 


1 


^(5) 


72 Q + 


33 D 


- 887 


56 


129 


8 


6x(3) - x^i) 


12 Q 


+- 6D 


- 26 


10 


17 


2 


6x(4) _ 2x^2) 


6Q H 


- 8D 


- 17 


13 


8 


1 


6x'-^^ - Sx^^^ 


68Q 4- 


30 D 


- 732 


52 


120 


9 
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With the nuUspace computation we now understand the constant / (Aq in [7 ). 
Thus for the minimal order ODE, one should have / = 1 since the minimal order 
ODE is unique. Setting Q = q and D = d + Dapp, where Dapp is the degree of the 
polynomial whose roots are apparent singularities, one obtains the exact relation 

Dapp = (d-l)(q-l) - C -1 (B.4) 

between the constant C and the degree Dapp. For x^^'^ one has d = 12, g = 7 and 
C = 37 giving Dapp = 28 which is [TUl [H] the degree of the polynomial carrying 
apparent singularities in the linear ODE of x^"^-*. For x^^-* one has c? = 7, g = 10 and 
C = 36 giving Dapp = 17, which is [TT] the degree of the polynomial carrying apparent 
singularities in the ODE of x^"^- Similarly for x^^\ with d = 72, q = 33 and C = 887 
we obtain the degree of the apparent polynomial Dapp = 1384, which is in agreement 
with what appears in the linear ODE for x*-^^ reduced to its minimal order. 

Note also that (|B.4p is valid for linear ODEs without an apparent polynomial 
(x'^-* and x^^'')- But there are cases where the parameter C is negative while the 
linear ODE has an apparent polynomial. This is the case we consider now. 

Appendix B.l. The ODE formula for the factors 

We first show how the apparent polynomials occur in a factorization of linear 
differential operators such as 

£ = LR 

where the factors L and R are monic0 and of minimal order, denoted respectively 
ql and q^. Denoting by Papp^ the apparent polynomial occurring in £, one knows 
that this polynomial should appear as an apparent polynomial in the left-operator L. 
It may happen that the right operator R also contains a polynomial Q of apparent 
singularities and this polynomial should not appear in C. For this to happen, the left 
operator L must have the roots of Q as true singularities. Furthermore, Q should 
occur in L to the power of the order of L, i.e. as Q^^ . The local exponents for L 
at any root of Q are — 1, 1, 2, •••,(7^ — 1. If we remove the singularity from L, 
the new linear differential operator L will have Q as an apparent polynomial and will 
occur as Q*^^^ with local exponents 0, 2, 3, • • • , qL- 

Consider, as an example, the series S for x^^"* annihilated by a seventh order 
linear ODE with L7 as the corresponding linear differential operator. We know that 
this operator factorizes as (among other factorizations p9)) ) 

Lt = L R = {Ah ■ Y^) ■ (Z2 • TVi). (B.5) 

Assume that the right-operator R is known. The aim is to produce the left-operator 
L by acting on S with R. The series R{S) will satisfy a linear ODE corresponding to 
L. 

For the right-operator R = Z2 ■ Ni, the left-hand side of the ODE formula (|B.3p 
reads 

dR-Q +qB- D ^Cr = 8Q+3D-9. (B.6) 

Putting these values into (|B.4p we obtain D^^^ — 4 which is the degree of the apparent 
polynomial Q occurring in R = Z2 ■ Ni . 



tl Normalization of the head polynomial of the linear differential operator. 
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The linear ODE for the left-operator L produced from the series R {S) when R 
is taken monic and of minimal order, has the ODE formula 

dL-Q +qL-D - Cl = 15Q+4D-1. (B.7) 



The degree of the apparent polynomial for L — Mi ■ Y3, computed by (jB.4p . is 
D^pp = 40, which is the degree of Papp (the apparent polynomial of L7, see above) 
plus three times the degree D^^p, and we still have the roots of Q appearing with 
multiplicity one in = 15. 

In computations modulo a prime, and for high order linear ODEs, it is obvious 
that it is easier to work with non-monic operators. This results in removing the pole 
part of the polynomial Q, leaving its apparent part in the left-operator L. 

As an example, we will reproduce the series R (S) with R non-monic but still 
of minimal order. The left-hand side of the ODE formula (|B.3[) . corresponding to L, 
reads: 

dL-Q +qL-D - Cl = 4Q + AD+32. (B.8) 

From (jB.4p we obtain the degree of the apparent singularities in L as 40 = 28 -f 3 x 4. 
Furthermore, recalling the left-hand side of the ODE formula (jB.31) (see Table [BT|) 



dQ + qD-C = 12Q + 7D-37 (B.9) 
for the full L7, one has d = du + dL and 

C = Cl +^— ^^-Cfl (B.IO) 

qB. - i 

+ ■ {{q - - 1) ^fpp + dqn -2qR + q-d 

Remark: Even if the various parameters in the ODE formula are now 
understood, we should recognize that we still do not know how this ODE formula can 
be proved, nor where it comes from. The various formulae dealing with the apparent 
polynomial degree (in fact upper bounds, e.g. [221 1301 131] ) in Fuchsian hnear ODEs 
introduce ingredients that go beyond our experimental mathematics framework. 

Appendix C. Some linear differential operators in exact arithmetic 

The linear differential operators V2, F2 and occurring in the decomposition of Ln 
ill = {Z2 ■ Ni) ®V2® (F3 • F2 ■ LI), (Cl) 
read respectively 

j3 + 8«^+W^ 1 + 7^ + 4^^ 

^ (l + 4w)(l-4w)u; (l-4w)(l + 4w)2u;2' ^ ' 



and 



F2 - D^ -§-D^ (C.3) 



witlfl: 

P2 = (1 -4w) - Pa, 



P2 ^ P2 



t Note that the factors (1 -|- 2ui) and (1 — ui) appear to the power one in both P2 and Pi. Linear 
differential operators can be Fuchsian without having descending powers of the factors giving rise to 
the singularities. 
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P2 = w • (1 - 4w) (1 + 4u)) (1 - w) (1 + 2w) (1 + 3-if + 4^2) 
(1 + w - 24 - 145 - 192 w* + 96 + 128 w^), 



Pi = (1 -w){l-4w){l + 2w)(^ 



40960 + 24576 + 51712 ~ 66816 



- 138176 - 88704 - 29940 - 5394 - 272 + 92 + 1 1 ^ + l 

Po = 262144 _ 65536 w^^ + 335872 - 934912 w^" - 743424 + 703488 w** 
+867776 w'^ + 371848 + 96744 + 14710 + 2144 + 398 + 9 - 11. 

It is possible to get rid of the apparent singularities occurring in F2, by multiplying 
F2 at the left, by a first order linear differential operator Ci, 

1 /go 256352914629 X 

73326885520'Vp2 l+3M;+4w;2y'' ^^■^> 

with: ^ ^ Aln(i?(u;)), 

where R{w) is a rational function (with integer coefficients) and 

go = 244820905584+ 1372135276587 w;+ 1384232623846 w^- 13621658367235 
-150856196156313^''- 1054439469518747 w'"^ - 3472090747016314 
-3873078043825712 + 3114022565962720 + 12058813946690432 
+10882841933451520 w^" - 1075293814167552 - 3544662480211968 w^^ 
-9348606615093248 

thus yielding a third order desingularized Fuchsian operator. This is the so-called 
"desingularization" procedure which preserves the Fuchsian character of the linear 
differential operators. Note however that the desingularization procedure does not 
preserve the remarkable property of global nilpotence of the highly restricted second 
order differential operator F2 . The new desingularized third order differential operator 
is no longer globally nilpotent because the first order differential operator Ci is 
not globally nilpotent. The breaking of global nilpotence comes from the factor 
256352914629/(1 + iw + 4^^) in (IC.4jl which is not a logarithmic derivative of a 
rational function. 

Next we focus on the "physical" singularity w — 1/4. One can change the 
operator F2 into a slightly simpler one as follows: 

/ \-ll/4 

F2 — > F2 ^ F2-[l- 4«;j . (C.5) 

where the dot corresponds to a multiplication of (differential) operators. It is 
important to note that the solutions of F2 around the "physical" singularity w — 1/4 
are in fact Puiseux series in u = [w — 1/4)^/^. In other words F2 rewritten in terms 
of the variable u is not singular aX w — 1/4. 

The calculations performed on Z2 yielded a modular form interpretation of Z2 
(see [B]). A crucial step corresponded to discovering the covering 

(l-w)(l-4w) ^ ' 

or Q(i, w) = 0, with Q{t, w) ^ At ■ - {bt - S) ■ w + t, 

which wraps the singularities of Z2 onto the three singularities of 2F1, namely 0, 1, 
00. Do note that the apparent polynomial for Z2 occurs as a vanishing condition of 
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the discriminant in w of the covering polynomial Q{t^ w): 

discrim(Q(t, = (9t-8)(t-8) = 8 ^ T, ^ , ^ - (C.7) 

Trying to perform a similar calculation for Fi in order to discover some modular form 
interpretation for _F2, we observe that it is not straightforward to find a covering, 
such as (jC.6p . wrapping all the singularities of F2 onto 0, 1 and 00, and such 
that the discriminant in w (like (|C.7[) ) of the corresponding covering polynomial 
Q{t, w), could correspond to the quite involved apparent polynomial of F2, namely 
1 + w - 24 - 145 u;^ - 192 w"^ + 96 + 128 w^. For these reasons we have failed 
in finding a modular form interpretation of the highly restricted linear differential 
operator F2. 

The third order linear differential operator F^ reads: 

= + {l + 2w) ^-Bl +2Ps^-D^ + ^, (C.8) 

where 

Ps = -w • (1 - 4'u;)(l + 4u;) (1 + ?« - 24u)^ - 145w^ - 192w^ + 96w;^ + 128w'^), 

P3 = {1 - w) {1 - 2w) {1 + 3w + 4w^) {1 + 2w)^ ■ P^ ■ p3, 

P3 = 5629499534213120^^^ + 5348024557502464^3^ - 62874472922742784^3^ 
+339080589913096192 ly^" + 132348214635397120 w^^ + 354600746294968320 w^^ 
+1383732497338073088^31 - 269118080922157056 u;^^ - 1021414905992970240 
+401943021895024640 + 378516473892569088^^7 _ 379126125978189824 
-181955521970962432 + 182991453503356928^^4 + 119809766351437824^^3 

-34528714733649920^2^ - 46719523456286720 - 1865897472688128 u;2° 
+9861412040736768^1^+ 1690374175916032^1*- 1285664678690816^1^ 
-304716171767808^1^ + 112170181177344^1^ + 30517814178816 w^'^ 
-7815766123264 wi3 - 2274047571904 w^^ + 456062896896 u;" 
+150282885872 u-i" - 10690267808 - 6048942832 - 486602112 
+33772908 + 25075632 + 4670454 + 13440 ^3 - 69066 - 5169 w - 63, 

P2 = 2582544170319337226240w^i + 2029141848108050677760 w'^" 
-32885932997405703143424 + 193641813610004500971520 
+20426022066743356162048 w"*^ + 288714242895676430090240 w^'^ 
+618280187651267892346880 u;"!^ - 648919373873770257186816 u;'i'^ 
-863129472633247214075904 ^"13 - 1021011939308518347112448 w^^ 
-220333306036159265112064 w"" + 1659564100832816225320960 
+588473220873831600619520^39 - 1065067759683713707802624 w' 



-9030793760523344150528 ^3^ + 805481511795301371871232 w" 



^,38 

,36 

-122169749668787845595136^35 ^ 629129357422714417053696^34 
-87120833646056343339008^33 + 304015333576904250753024^32 
+143209349380404068483072 i(;3i - 67135556652765458464768 u;30 
-68161001548708224958464^29 - 1506006178531414900736 
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+15819782847593648750592 + 4086678104179764363264 w^e 

- 1909688698451711754240 tt;^^ - 970204468920909561856 
+94919087350092267520 w;23 + 123918740818141650944 

+4687965654930399232 - 10707547611722045440 
-1144659629046790144 + 806082415949659264 
+149774462467091328^1^ - 48096268859594016 w^*^ 
- 16578560990131776 w^^ + 424243043096032 w^'^ + 905149437225280 w^^ 
+139111711404072^1^ - 7972709043232 «;ii - 6405062530332 u;i° 
-1037367028148 + 6971928216 + 30288912150 w'^ 
+3873954375 - 115755798 - 60227304 
-3099678 + 211068 w'^ + 21432 w + 315, 

Pi = 19267255250108152471879680 + 26483031235932970358407168 
-256308802040991428795957248 w^^ + 1579949167665869307621933056 w^^ 
+650374789771441405855531008 + 5216643706804247528946532352 
+7917834014591751323461353472 w^^ - 3351835287019392824172871680 w^^ 
- 1 1064588131657140234556014592 w^^ - 34985695129493606924629835776 w^"^ 
-27150264881506217380601135104 ty^i + 21916196537425570804428439552 
+42001979686686526227299172352 w^^ + 25840385187494624677295292416 w^® 
-5424492229252674644950908928 w''^- 17794118224994570424773771264 w'"^ 
+2915867386820035753799581696 w^^ + 7186426242807565546487283712 w^^ 
-14774359259974734620101967872 w^^ - 14907706789958430400446464000 
+8071574290338795697467293696 w^^ + 15504536229837797153841348608 
+2626569595260883926907879424 w^^ - 6860129397955391680596148224 w^^ 
-4292305256193524038115524608^^^ + 797725481517011621914869760^3'' 
+1720922858808986948924866560 w^^ + 390533143866840910481326080 w^^ 
-290442221202989562927775744 w^^ - 165938014975046925940686848 w^^ 
+9585161397342427263533056 w;3i + 28502663270123757533921280 u;^^ 
+4489386799471924718338048 u;^^ - 2717486608897256297267200 
-908343774367384075960320 w^'' + 162310791240979996000256 w"^^ 
+103239269328878845726720^25 _ 7750369303138783333376^24 

-11489552784013679223808 w;23 _ 712249616867767788544 
+991748945187237072640^21 + 252693598182584513344 
-21130273995450588928^1^ - 20284808101979844832 
-3056348368556274592 u;!^ + 345270164930943040 

+205893879174875432 w^^ + 28654368006663856 

-2030520435693824 to" - 1374304588556840^12- 166988492206488 u;" 
+12760292849076 + 5484990319472 + 367504601004 
-50197207920 - 9218315844 - 277909095 
+48467763 + 5648070 + 265293 + 4620 w - 63, 
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Po = 69634127209802640463075737600w™+ 102981 137018052571618170896384 
-1033960403593443123509300559872wSS+ 6716228494346939277472100777984 
+830768383072984903026797969408 w^^ + 34119483032722461380174390755328 w^^ 
+37151403895216351475147854577664 w^^ + 7596402077224314128199487324160 w^^ 
-5774876585209674171391426953216010^2- 30949330267349771483063051 1968256 
-232460008226528101141464649564160w^°+ 23931702098177545910680337514496 
+427559960442089709631493273288704 w^^ + 767305599958046596665201651613696 
+238126263520324803598765665550336 w^^ - 489368606355635167460530948407296 w^^ 
-411394912009392610164715657625600 w^^ - 9190571673284172503536766025728 w^^ 
+4354675264071893610129679974400 w^^ - 197517596854575225398680040243200 w^^ 
- 100108534915833684237428566523904 + 237797067305428725186438474235904 w'^^ 
+297736260249409824018326167224320 w^^ + 2827068864630764422668662865920 w"*^ 
-212811801685085801466392370741248 w^'^ - 132552668920641958581606566330368 w^^ 
+36227143968974260317176152981504 w'^^ + 80360530046171440422025236054016 w'^^ 
+25044389743118121276003435151360 w^^ - 16775611588713607092922243612672 w^^ 
-14375299221388934261580907937792u;''°- 797769185002542420001267122176 
+3060062577366019941153762181120w^^+ 1048552246961478552732246736896 
-286300750377610217893186764800 w^^ - 238703363798670426041267257344 w^^ 

,33 

+6256946928524452096094240768 w^"^ - 3100756305550863462745636864 w^^ 
-1271279614733459684395712512 + 173570224057030875371798528 
+187963836513544173604265984^28 + 20098454205158749911726080^27 
-15220510128449109866076160 _ 5530208120756891132317696^25 
-31540092813892142535680 w^^ + 410046209080624124809344 w^^ 



-6453603219285212538454671360 w^* + 32067040600375745464846254080 w 



+95695972411021353163264^22 _ 3799963752408310388096^21 
-6052638686215258044992^20 - 1044163538474290733536 
+69113265719111269072^1® + 55225911443186243360 u;^^ 
+6988584609018020432 uii^ - 714608406420145560 u;^^ 
-313788688846958472 wi'' - 23383932527942400 
+4392065243452176^^2 ^. 942992856333120 + 18782060660376 
-11352161581890 - 1093090772088 w® + 23284774974 zz;^ 
+9267369222 + 542276796 - 59916 
-3757362 - 465618 w'^ - 20622 w - 126. 
We note that F3 can be simplified as follows: 

F3 — > F3 ^ Fs--, where: (C.9) 

^ = u;2 . (1 - 4u;)^/2 . ^ 4^^/^ • (1 - w) • (1 + 2w) 

X (l + 3w + 4^2) . App(F2), (C.IO) 
where App{F2) denotes the apparent polynomial for F2, namely 1 + w — 24w2 — 
145 w'^ — 192 w'^ + 96 + 128 w'^ , and where the dot in (|C.9p denotes the multiphcation 
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of (differential) operators. This just amounts to multiplying the solutions of by /z. 
Remarkably F^ is no longer singular aA, w = I nor at the two roots of the quadratic 
1 + 3w + Aw'^ — 0. We find the following exponents at the remaining singularities: 

w = 0, 0, 1, 3, (log2), 

w = 1/4, 0, 1, 3/2, 

w = -1/4, 0, 1, 5/2, 

w = oo, - 18, -18, -16, (log^), 

w = 1/2, 0, 1, 1/2, 

w = -1/2, 0, 1, 1/2, 

App{F2) = 0, 0, 2, 3. 

We have here an illustration of what we described in [Appendix B[ l where the third 
order linear differential operator F3 reads: 

F3 = Ps,ng ■ App[F2f ■ App{Fs) -Bl + ■■■ (C.ll) 

where Psing denote the "true" singularity polynomial of F3. We remark that the 
apparent polynomial of F3 is the apparent polynomial appearing in the product 
F5 — Fy, ■ F2- The polynomial App{F2) is the apparent polynomial of F2- It appears 
at the power of the order of F3 for which it is a pole. When rescaled as done in F3 
the roots of App[F2) become apparent singularities of F^. 

Note that the formal series of the linear differential operator ^3 are Puiseux series 
around all the singularities except w — Q and w = 00. These are the only singular 
points around which F3 has logarithmic solutions. When the third order operator F^ is 
rewritten in terms of the variable u = {w — WsY^'^ , where Ws is any singularity other 
than w = OT w — 00, F^ is no longer singular at (in particular the ferromagnetic 
critical point w — 1/4 is no longer singular in the variable u — {w — 1/4)^/^). 



Appendix D. Experiment: rational reconstruction of the apparent 
polynomial in F3 

Write the linear differential operator F3 as 

F3 - V3iw)Pappiw^')-Bl, +p2iw)P2{w'')-Bl 

+ Viiw)Piiw'^')-Bu, +Poiw^") 

where Vi{w) account foiQ the known multiplicities, and the argument in the 
polynomials is used to show their respective degrees n. Assume that this linear 
ODE has been obtained for many primes. We want to carry out the rational 
reconstruction for each polynomial separately, basically because the polynomials at 
the lower derivatives are harder to obtain. 

As it comes from our solver, the polynomial Papp cannot be reconstructed with 
nine primes. If we multiply all the mod prime coefhcients by 2"^^ the rational 
reconstruction will be successful with eight primes. If we multiply by 2^"^ the 
reconstruction succeeds with six primes. It should be noted that when the number of 
primes is not sufhcient, the correctly reconstructed coefficients will be those of lower 
degrees or higher degrees depending on the magnitude of the scale used to multiply 
the coefficients. This then calls for a scaling of the variable itself. If we change the 

X These Pi{ui)'s are different from the ones in l|62|l . 
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variable w to w/2 and multiply all coefficients by 2^°, the rational reconstruction is 
successful with just five primes. It is fortunate that the apparent polynomial is the 
easier polynomial to reconstruct. It will be used in further checks. 

How can one guess the scaling (e.g. 2^^ and 2^") mentioned above? We have 
found that 2^^ is the magnitude of the lower coefficient in Vsiw), which is an exactly 
known polynomial. The scaling 2*° is around the magnitude of the lower coefficient 
in 7^3 (w) • Papp- More than an educated guess, we have an almost deterministic 
procedure to find the proper scaling factors to improve our rational reconstructions. 
This experiment shows that the rational reconstruction is actually easier when the 
underlying physical problem is taken into account, leading to proper scaling factors. 
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